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Abstract. Given a pair of distinct vertices u,v in a. graph G, we say 
that s is a junction of u, v if there are in G internally vertex disjoint 
directed paths from s to it and from s to v. We show how to characterize 
junctions in directed acyclic graphs. We also consider the two problems 
in the following and derive efficient algorithms to solve them. Given a 
directed acyclic graph G and a vertex s £ G, how can we find all pairs 
of vertices of G such that s is a junction of them? And given a directed 
acyclic graph G and k pairs of vertices of G, how can we preprocess G 
such that all junctions of k given pairs of vertices could be listed quickly? 
All junctions of k pairs problem arises in an application in Anthropology 
and we apply our algorithm to find such junctions on kinship networks 
of some brazilian Indian ethnic groups. 



1 Introduction 

Given a directed graph G, a vertex s € G is a junction of a pair u,v £ G, 
u ^ V, if there exist two internally disjoint directed paths from s to u and 
from s to V. A vertex s G G is a lowest common ancestor, LCA for short, 
of a pair of vertices u,v, u ^ v, if it is a junction and for every junction 
s' s of the pair u,v, there is no directed path in G from s to s'. In rooted 
trees, junctions and LCAs are the same vertices and given a pair u,v £ G, 
u ^ V, the LCA is unique. On the other hand, in directed graphs, acyclic 
or not, we can have several different LCAs or junctions for a certain pair 
of vertices. The problem of finding out whether s is a junction, or an LCA 
of a pair n, v can be done efficiently using a maximum-flow algorithm. So, 
we can determine whether a given vertex is a junction, or an LCA, on 
directed graphs in polynomial time. There are faster algorithms that find 
out junctions and LCAs in rooted trees and directed acyclic graphs. 
Related Works. A problem that have been studied in rooted trees and 
directed acyclic graphs (DAGs) is how could we preprocess a given graph 
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such that a query to a representative LCA for any pair of vertices could 
be done quickly. We have found in the literature fast algorithms that 
solve this problem. Aho, Hopcroft and Ullman (I1973P have shown how to 
preprocess rooted trees with n vertices in 0{na{n)) time, where a{n) is 
the number of times that we must apply log2 to n to obtain a number less 
than or equal to zero. Queries are answered in constant time. Harel and 
Tarjan (|1984p have shown how to preprocess in 0(n) time and to answer 
queries on the representative LCA in constant time. Other algorithms to 
preprocess rooted trees appeared in the works by Berkman and Vishkin 
([EMD, Nykanen and Ukkonen (fTOMll and Wen (|T9M1) . The equivalent 
problem to DAGs is known as ALL-PAIRS-LCA (Bender, Farach-Colton, 
Pemmasani, Skiena and Sumazin 120050 . Given a DAG with n vertices, 
the preprocessing phase proposed by Bender, Farach-Colton, Pemmasani, 
Skiena and Sumazin (pnn5]) spends d{n^-^^^Q and a query to a represen- 
tative LCA can be answered in constant time. Later, Czumaj, Kowaluk 
and Lingas (12007P have solved this problem to DAGs with preprocessing 
time 0{'n?-^'^^) and Eckhardt, Miihling and Nowak (|2007p with expected 
time 0{n^ log n). 

In the problem ALL-PAIRS-ALL-LCAS we want to know how to pre- 
process a given DAG with n vertices and m arcs such that a query 
to all LCAs of a pair of vertices can be done quickly. Baumgart, Eck- 
hardt, Griebsch, Kosub and Nowak (j2007p have developed algorithms 
that solve ALL-PAIRS-ALL-LCAS. The upper bound of the preprocessing is 
0(min{n^m, n^'^^^}). Eckhardt, Miihling and Nowak ()2007p have devel- 
oped two algorithms for this problem, one with preprocessing expected 
time 0(n^ log log n) and the other with preprocessing time 0(n^'^^^^). 

Yuster ()2008p considers in his work the all-pairs-junction problem, 
i.e., how to preprocess a DAG with n vertices such that a query to a 
representative junction can be done quickly? Any algorithm that solves 
ALL-PAIRS-LCA could be used to solve all-pairs- JUNCTION but Yuster 
()2008j) have shown that we can preprocess in 0{n^) time and answer 
queries in constant time, where w < 2.376 is the exponent of fast boolean 
matrix multiplication. 

In the problem treated by Tholey (j2005p there is given a DAG G, 
si,S2 £ G and k pairs of vertices {ui,vi), . . . , {uk,Vk)- For each pair 
{ui,Vi), I < i < k, a tuple (si, ti, S2, ^2) with {ti,i2} = {ui,Vi} is 
printed out in constant time if there are two disjoint directed paths 
from si to ti and S2 to ^2- Before that he constructs a modified ver- 



/(n) = 0{g{n)) if there is a constant c such that f{n) — 0{g{n) log"^ n). 



3 



sion of the data structure proposed by Suurballe and Tarjan ()1984p in 
0(nlog2 n + {m + k) log2+^^+k) / (n+k) n) time. 

Our Application. This work has been motivated by an apphcation that 
arises on the field of Anthropology. It is given a kinship network (a DAG) 
that models the parent-child relationships from a determined society (in 
our case Indian Brazil ethnic groups). Moreover, the set of weddings 
among individuals from that society is given. We want to know if the 
partners of the wedding are relatives in some degree. From the anthropo- 
logical point of view, the junctions, i.e., common ancestors with disjoint 
descendant lines, are the objects of interest (dal Poz and Silva l2009|) . 
In this case, we want to find all junctions of many pairs of individuals of 
the network. Thus, we treated here the problem /c-pairs-all-junctions, 
that is, given a DAG G and k pairs of vertices of G, how can we prepro- 
cess G such that a query to all junctions of k given pairs of vertices can 
be done quickly? 

Main Results. We develop an algorithm to solve the following problem 
named single- junction- all-pairs. Given a DAG G with n vertices 
and m arcs, and a vertex s € G, construct a data structure that allows 
to find all pairs of vertices for which s is a junction. The time spent by 
it is 0(m). We can use another data structure named dominator trees 
(Aho and Ullman, I1972p to solve single-junction-all-pairs problem. 
We construct a dominator tree T rooted in s in linear time considering 
the graph induced by the vertices descendants of s in G. The pairs u, v 
that have s as a lowest common ancestor in T are the pairs that have s 
as a junction in G. However, our data structures is simpler than the data 
structures for dynamic LCA queries in trees (Cole and Hariharan, I2005P 
used to construct a dominator tree in DAGs. 

To solve the k-pairs-all- junctions problem, we just solve the prob- 
lem single- junction-all-pairs for all s in G listing (or storing) s to 
pair u, V, if s is a junction of n, v. The time spend by it is 0{n{m + k)). 

The algorithm to K-PAIRS-ALL- JUNCTIONS problem can be used to 
solve the k-pairs-all-lcas problem. Given a DAG G with n vertices 
and m arcs, and k pairs of vertices, we first find the transitive closure of 
G, find all junctions of the k given pairs of vertices, and finally for each 
given pair u, v and for each junction Sj G Jiu, v), we verify if there exists 
an arc Si — sj in the transitive closure of G, where Sj ^ Si and Sj € J(m, v). 
If yes, then Si cannot be an LCA of the pair ti, v because Sj has Sj as a 
common descendant to v. If no, then Sj is an LCA of the pair u, v. 
Thus, the time spent by this method is 0{n^ -|- n(ra -|- /c) -|- n^/c), where 
O(n^) is the time of transitive closure of G. If A; = o(n^'^^^^), then our 
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simple approach for K-PAIRS-ALL-LCAS is faster than the best algorithm 
(worst case) known for solving ALL-PAIRS- ALL-LCAS problem. 

Any algorithm that aims to list all junctions of all pairs of vertices 
spends !7(n^) time, in the worst case. To see that consider the example 
in Fig. [H where each vertex in the first line is a junction of every pair 
of vertices in the second one. Thus, to list all the junctions of all i7(n^) 
pairs would spend I2(n^) time. 

Q{n) 
n{n) 

Fig. 1. Pairs witii i7(n) junctions. 

Organization. This text is organized as follows. In Section 2, we intro- 
duce some notation used in this work, recall an ordinary search algorithm 
(depth-first search) and discuss some simple data structures. In Section 
3, we characterize junctions in DAGs. In Section 4, we describe the al- 
gorithm for SINGLE- JUNCTION-ALL-PAIRS problem, and in Section 5, we 
make our last considerations. 

2 Initial Concepts 

Given a DAG G and vertices s,t G G we say s is a parent of t or t is 
a child of s, when there exists an arc s — t € G. We denote the set of 
parents of a vertex u in G by 6q (u) and the set of children of a vertex u 
in G by 5q{u). We say s is an ancestor of t or t is a descendant of s when 
there exists in G a directed path from s to t. Usually, we treat a directed 
path P = {s = tQ,ti, . . . , tk-i = t} as an ordered set of vertices with the 
property that U — tj+i G G, for alH = 0, . . . , A; — 2. 

If a vertex s is a junction of the pair of distinct vertices u,v, then 
there are directed paths internally disjoint P from s to u and Q from s 
to V. Note that, if s is an LCA of u, v, then s is a junction of u, v, but the 
converse is not always true. As we mentioned before, a pair of vertices 
u,v G G can have many LCAs and junctions. We write J{u, v) to denote 
the set of vertices of G that are junctions of the pair u, v. The sets J(n, v) 
and J{v, u) are equal. If u is a proper ancestor of v, then u € J{u, v). We 
suppose that J{u, u) = 0, for all u G G. Analogously, LCA{u, v) denotes 
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the set of LCAs of the pair u, v. If we want to specify the graph G we 
are working on, then we write LC Aq{u,v). When G is a rooted tree, we 
abuse of the language writing s = LCAg{u,v). 

We call an arborescence a rooted directed tree, such that: 1 - there ex- 
ists a unique vertex with in-degree equal to zero called root; 2 - all vertices 
except the root have in-degree one; and 3 - there exists a directed path 
from the root to all vertices in the arborescence. An arborescence with 
root s is denoted by Tg. Any vertex n G is the root of a subarborescence 
denoted by Tu ^Tg. 

Given a DAG G and a vertex s € G, we use a depth-first search to 
construct an arborecence Tg composed by root s and all descendants of s 
in G. We denote by Fg the set of all arcs of G with both ends in Tg that 
are not arcs of Tg. For each vertex u & Tg, we maintain the integer post- 
order value stored in array post indexed by vertices. We also maintain the 
integer values minpost[u], that is, the minimum post-order values in Tu- 
After such construction, we can make a partition of the arcs with both 
ends in Tgi 1. Arcs of the arborescence Tg] 2. External descendant arcs 
- are the arcs in Fg of the form s — v, v G Tg. Note that v is a child 
of s in G but is not a child in Tg; 3. Internal descendant arcs - are the 
arcs u — v(zFg,u^s, such that there exists a directed path in Tg from 
u to v; 4. External crossing arcs - are the arcs u — v G Fg such that 
u G and v G T^^, where si and S2 are two different children of s in 
Tg-, 5. Internal crossing arcs ~ are the arcs u — v G Fg such that u G Tg-^, 
V G Tg-^ , si is a child of s in Tg and there are no directed paths in Tg from 
u to V neither from v to u. A similar partition of the arcs has been done 
by Sedgewick (I200ip and Dasgupta, Papadimitriou and Vazirani (I2006p . 
Such construction has some properties. 

Property 1. For any arc u — v with both ends in Tg we have post[u] > 
post[v]; If si and S2 are two children of s in Tg and post[si] < post[s2] 
then, there is no directed path from p to q, where p G Tg-^ and q G Tg^; 
and post[p] < post[q], for all p G Tg^ and for all q G Tg^- 



Proposition 1. Let G be a DAG, s a vertex in G and Tg an arborescence 
constructed by depth-first search as before described. Let si be a child 
of s in Tg and u be a vertex in Tg-^. Let P = {s = uq,ui, . . . ,Uk-i = 
u} be a directed path from s to u in G. If the vertex Ui, for some i = 
1, . . . , A; — 1, is the first vertex of P in Tg-^, then all remaining vertices of 
P, Ui+i,. .. ,Uk^i, belong to Tg^. □ 
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3 Some Results for Junctions in DAGs 

In this section, let us consider a DAG G, a vertex s E G and an arbores- 
cence Tg constructed as mentioned before. The arcs of G that do not 
belong to Tg or Fg are not used in any directed path beginning in s. So, 
they are not important to obtain the pairs u, v that have s as a junction. 
The following proposition can be quickly checked. 

Proposition 2. Let Tg be an arborescence with root in s, si and S2 dif- 
ferent children of s in Tg and two vertices u,v € Tg. If u G Tg^ and 
V G Ts2, then s G J{u,v) and s = LGAt^{u,v). □ 

It remains to us to find out whether s is a junction of pairs u, v when u 
and V belong to a same subarborescence. Note that if we take two vertices 
u and V that belong to subarborescence Tg. , s j child of s in , and if there 
is no external arc p — q with q G Tg., then s cannot be a junction of this 
pair, since all directed path from s to n and from s to u necessarily share 
the vertex Sj. Therefore, if s G J{u,v) with u,v (z Tg., then there exists a 
pair of internally disjoint paths, from s to n and from s to v, and one of 
them has an external arc. 

Additionally, we note that for any vertex si child of s in Tg and 
a pair u,v in T^j, if s G J{u,v), then there exists a pair of internally 
disjoint paths in G from s to n and from s to v, such that all internal 
vertices from one of them are in T^^ . We can prove this by contradiction. 
Consider two internally disjoint paths in G, P = {s = uq, . . . , u^-i = u} 
and Q = {s = vq, . . . ,vi-i = v}. Suppose that P and Q do not have 
all vertices in Tg^. Consider the first vertices Ui and Vj of P and Q, 
respectively, that belong to Tg^ . By Proposition [H the directed paths 
Pi = {ui,Ui-^-l, . . . ,Uk-i} and Qi = {vj,Vj^i, . . . ,vi-i} are internal to 
Tg^. Take the directed path P2 from si to Uj in T^^. If P2 H Qi = 0, then 
we can construct a directed path from s to u ({s} U P2 U Pi) with all 
vertices in , where {{s} U P2 U Pi) n Q = {s}. So, P2 n Qi / 0. Cah v' 
the vertex from this intersection nearest to vi-i. In the same way, take 
the directed path Q2 from si to Vj in T^^. If (52 H Pi = 0, then we can 
construct a directed path from s to f ({s} U Q2U Qi) with all vertices in 

and {{s}UQ2UQi)nP = {s}. So, Q2nPi / 0. Cah u' the vertex from 
this intersection nearest to Uk-i- Thus, we find a cycle with the directed 
paths from Vj to v' , from v' to Ui, from Uj to u' and from u' to vj. A 
contradiction, since G is a DAG. Therefore, all internal vertices to P or 
Q are in T^^ . This fact helps us to prove the following lemma. 
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Lemma 1. Let si be a child of s in Tg, u and v belong to Tg^ and z = 
LCAt^ {u, v) with z u,v. The vertex s € J(n, v) if, and only if, z belongs 
to a pair of internally disjoint paths from s to u or from s to v. 

Proof. Take a pair of internally disjoint paths P = {s = uq, . . . , ttfc-i = u} 
and Q = {s = vq, . . . ,vi-i = v}. Suppose all internal vertices of P are 
in Tg-^ and Q has an external arc Vi — Uj+i entering in T^^. Consider 
R = {si = zq, . . . , Zm-i = z} to be the directed path in Tg^ from si to 
z. Consider P' = {z = Uq, . . . , u'^_i = u} the directed path in T^^ from z 
to u and Q' = {z = v'q, . . . , v'^_i = v} the directed path in Tg^ from z to 
V. Here is really important to note that, by the construction of T^, any 
path that passes through a vertex in Q' {P'), and then passes through a 
vertex in P' {Q'), cannot come back to Q' {P'). Let us divide the proof 
in two cases. 

Case 1. There is no vertex of Q in R. If there is no vertex of Q in P', then 
we can make Z = {s} U RU P' and we have Z and Q internally disjoint. 
If there is some vertex of Q in P', then consider Vj = u[ the vertex in 
Q n P' nearest to z. So, the directed paths U = {s = vq, . . . , Vj} U {vj = 
u'^, . . . u'f^_i = u} and Z = {s} URUQ' are internally disjoint. See the two 
first illustrations in Fig. [2j 

Case 2. There are vertices of Q in R. Let us consider that the vertex 
nearest to z in (P U Q) Ci R belongs to Q, i.e., Vw = zt- The case when 
this vertex belongs to P is symmetric. 

Case 2.1. There is no vertex of P in Q' . Then, the following directed 
paths are internally disjoint: P and Z = {s = vq, . . . ,Vw} U {vw = 

Zt, . . .,Zm-l = Z}U Q'. 

Case 2.2. There are some vertices of P in Q'. Consider Up = v'^ the 
vertex in P n Q' nearest to z. So, the following directed paths are inter- 
nally disjoint: Z = {s = vq, . . . , v^} U {v^ = zt, ■ ■ ■ , Zm-i = z} U P' and 
U = {s = Uq, . . . ,Up} U {Up = v'g, . . . , v'^_i = v}. 

Therefore, if s € J{u,v), then there exists a pair of internally disjoint 
paths with z belonging to one of them. This ends the first part of the 
proof. 

The converse is easier. By definition, s is a junction of the pair n, v if 
there is a pair of internally disjoint paths, independently if z belongs to 
one of them. □ 

The Lemma[2]is a consequence of Lemma[TJ It characterizes s as a junction 
of a pair u,v when we have z = LCAx^{u,v) and z ^ u,v. Lemma [3] 
characterizes s as a junction of a pair u,v in the case where z = u oi 
z = V. 
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Fig. 2. New internally disjoint paths Lemma [T]- Cases 1 and 2. One of them includes 
vertex z = LCAts{u,v). The vertex £ Q is the nearest to z in {P U Q) D R. The 
unlabeled vertices are not important. They are in the figure just to indicate the input 
vertices in a path of T,. 

Lemma 2. Let si be a child of s in Tg, u and v belong to T^j, z = 
LCATg{u,v) and z ^ u,v. The vertex s € J{u,v) if, and only if s € 
J{z,u) or s ^ J{z,v). □ 

Lemma 3. Let si be a child of s in Tg, z and u belong to Tg^ and z he 
a proper ancestor of u in Tg, i.e., z = LCAt^{z,u), z ^ u. The vertex 
s € J{z,u) if, and only if, s € J{z,t), for some t € 5q{u). 

Proof. To prove the first part we take two internally disjoint paths P = 
{s = no, • • • , Uk~i = u} and Q = {s = zq, . . . , zi^i = z}. We know that 
Uk-2 € <^g(^) since it is a parent of u in G. Therefore, the directed paths 
U = {s = uq, . . . , Uk-2} and Q are disjoint, and then s G J{z, Uk^2)- 
To prove the converse, we consider two internally disjoint paths P = 
{s = to, • • • , ifc-i = t} and Q = {s = zq, . . . , zi^i = z}. If u = tj, for any 
i = 1, . . . , k—2, then the directed path {u = tj, tj+i, . . . , tk-.i,u} is a cycle. 
But G is a DAG. Thus, u ^ P . li u = zj , for any j = 1, ... ,1 — 2, then u is 
a proper ancestor of z. By assumption, z is a proper ancestor of u. Thus, 
we can again produce a cycle contradicting the fact that G is a DAG. 
Then, u ^ Q. Therefore, we can use the arc t — u extending the directed 
path P and constructing two internally disjoint paths U = P L) {u} and 
Q. So, s S J{z, u). □ 

The next two lemmas help us to detect, respectively, some pairs of vertices 
that have and some that do not have s as a junction. 

Lemma 4. Let si be a child of s in Tg, z and w belong to Tg^ and z be 
a proper ancestor of w in Tg, i.e., z = LCAt^{z,uj), z ^ w. Consider 
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the directed path from z to w in Ts, Z = {z = w'q, . . . , w'f^_^ = w}. If the 
vertex s ^ J{z, w^), for all i = 1, . . . , k — 2 and s € J{z, w), then 

a. the vertex s G J{w'^, w), for all i = 1, . . . ,k — 2; and 

b. for all pairs of vertices u,v, u (z Tz \ and v G T^, we have s G 
J(u, v). 

Proof, a. Consider two internally disjoint paths P = {s = wq, . . . , wi-i = 
w} and Q = {s = zq, . . . , Zm-i = z}. We know that {Q (1 Z) \ {z} = 0, 
since Z\{z} contains proper descendants of z and Q\{z} contains proper 
ancestors of z. Note also that (PCiZ) \ {w} = 0. If it is not the case, then 
s G J{z,w'^) for some i = — 2, against our assumption. So, we 

can extend the directed path Q using the path Z from z = w'q until w'^, 
for all i = 1, . . . , A: — 2. Therefore, we can construct new directed paths 
Ri = Q U {z = w'q, . . . , w'j}, for alH = 1, . . . , A; — 2, such that the paths 
Ri and P are internally disjoint and s G J{w'j^,w), for alH = 1, . . . , A: — 2. 
b. Take u G Tz\Tw and v G Tw Let w'j = LCAt^{u,v). By Lemma H] a, 
we have, in particular, s G J{w'j,w). Therefore, there exist directed paths 
internally disjoint P from s to w'j and Q from s to w. Note that we can 
extend the directed path Q adding the path in Trfju from w to v and this 
new path, denoted by R, does not intersect P or Q\ {w}. If it does, there 
would be a vertex ancestor and descendant of w. Thus, R H P = {s}. 
Therefore, we have s G J{w'j,v). If it = w'j, then is done, li u ^ w'j, then 
by Lemma m s G J{u,v). □ 

Lemma 5. Let si he a child of s in Tg, z and w belong to T^^ and z be 
a proper ancestor of w in Ts, i.e., z = LCAt^{z,w), z ^ w. Consider 
the directed path from z to w in Tg, Z = {z = w'q, . . . , w'f^_^ = w}. If the 
vertex s ^ J{z,w'^) for all i = 1, . . . , k — 1, then s ^ J{w'^,w'j), for all 
i,j = l,...,k-l. 



Proof. Proof by contradiction. Suppose that for some i,j = l,...,k — l, 
we have s G J{w'^,w'j). Consider P = {s = xq, . . . ,Xk-i = w'^} and Q = 
{s = yo) • • • 1 Vi-i = w'j}. Consider the directed path on arborescence Tg, 
R = {s = z'q, . . . , z'^_i = z}. Define W as the set of vertices of i?n (PUQ). 
Take the vertex z'^ € W nearest from z. Suppose that z'^ = Ur G Q (the 
other case, G P is analogous). Then, the directed paths Z = {s = 
Uo, . . . ,yr = z'g} U {z'g, -Zg+i, . . . , = z} and P are internally disjoint. 
Therefore, s G J{z, w'^), a contradiction. Thus, for all i, j = 1, . . . , A: — 1, 
we have s ^ J{w'^,w'j). □ 
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4 Algorithms for Junctions in DAGs 

Let us develop an algorithm that solves the problem where it is given a 
DAG G and a vertex s ^ G and we want to know which pairs of ver- 
tices u,v & G, u ^ V have s as a junction. We use a data structure for 
disjoint sets represented by rooted trees and we apply the path compres- 
sion heuristic. In this data structure, each set is identified by an element 
called representative (chosen as the root of a tree). See Gormen, Leiser- 
son, Rivest and Stein (|200ip . for more details. We use an array p indexed 
by vertices to represent such structure. For each vertex v, p[v] points to 
its respective parent in the data structure for disjoint sets. If = v, 
then f is a representative. The Single-Junction- All-Pairs algorithm 
receives two vertices s,z G Ts and constructs a data structure for disjoint 
sets, represented by the global array p, such that s is a junction of a pair 
u,v E Tz ii, and only if, p[u] ^ p[v]. 

Single- Junction- All-Pairs(s, z) 
1 if \Tz\ > 1 



2 then w ^ vertex {post[z] - 1) 

3 while post[w] > minpost[z] 

4 do if 3 t G S(j{w) with p[z] ^ p[t] 

5 then plw] ^ w 

6 Single- Junction- All-Pairs(s, tL)) 

7 w -(r- VERTEX{minpost[w] — 1) 

8 else p[w] z 

9 w vertex (post [w] - 1) 



The function vertex(p) receives an integer p and returns a vertex 
V such that post[v] = p, if post[v] is defined; or a dummy vertex with 
post-order value equal to —1, otherwise. It spends constant time. 

To construct a data structure for disjoint set such that s is a junc- 
tion of a pair u,v (z Tg if, and only if, p[u] ^ p['v], we have to make the 
calls Single- Junction- All-Pairs(s, Si), for all Si G S^^{s). We assume 
that Tg has been constructed as before described and the array p has been 
properly initialized as follows. For all arcs u—v € T^, = n, p\s\ = s and 
p[si] = Si, for all Si G 6^^ (s). Initially we have (s)| -|- 1 disjoint sets and 
for now u, v with u G T^, and v G Tg^, Si and Sj different children of s in 
Tg are pairs of vertices for which s is a junction, as we have said in Propo- 
sition [21 Let us see that the algorithm Single- Junction- All- Pairs is 
correct. 
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Lemma 6. Suppose a call to Single-Junction- All-Pairs(s, Si), where 
Si € 6^^ (s) . In line 3, the following assertions are true for the parameter 
vertex z and for all pairs of vertices u,v G T^., such that post[u\,post[v] > 
post[w] + 1. 

— The vertex z is a representative; 

— The values p[u] and p[v] point to a representative; 

— For all vertices x £ Ts-\ T^, p\z\ ^ p{x\; and 

— The vertex s G J{u,v) if, and only if, p[u] ^ p\'o\. □ 

Note that each arc in G is explored once by Single- Junction- All- Pairs 
algorithm. So, the next result follows. 

Theorem 1. Given a DAG G with n vertices and m arcs, and a vertex 
s £ G, we can construct a data structure for disjoint sets in time 0{m), 
and answer in constant time for any pair of vertices u,v €z G, if s is or 
is not a junction ofu,v. Moreover, after of such construction, all k pairs 
of vertices that have s as a junction can be printed out in 0{k) time. □ 

Theorem [T] solves the SINGLE- junction- all-pairs problem and helps us 
to solve other problem associated. 

Corollary 1. Given a DAG G with n vertices and m arcs, a vertex s G 
G, and k pairs of vertices of G, {ui,vi), . . . , {uk,Vk), we can answer in 
0{m + k) time which pairs {ui,Vi), 1 < i < k, have s as a junction. □ 

Applying Corollary [J for each s e G, we solve the /c-pairs-all- junctions 
in 0{n{m + k)) time. We spend, respectively, 0{m) or 0{nk + m) space, if 
we print out during the execution or we store all junctions of the k given 
pairs. 

5 Conclusion 

Given a DAG G with n vertices and m arcs, and k pairs of vertices of G, 
we have interest in listing all junctions of the given pairs. This problem, 
called A;-PAIRS- ALL-JUNCTIONS, arises in some application in Anthropol- 
ogy. It is solved in two parts. First, we construct a data structure for 
disjoint sets in time 0{m) such that a pair of vertices u,v in G has s as a 
junction if, and only if, u and v are in different sets. So, from k given pairs 
of vertices, we find out which of them have s as a junction in 0{k) time 
since we can say, in constant time, if s is or is not a junction of any pair 
of G. Second, we apply this idea for all s in G. Therefore, to solve the k- 
PAIRS-ALL- JUNCTIONS problem we spend 0{n{m + k)) time and, or 0{m) 
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space, if the junctions are listed during execution, or 0{nk + m) space, in 
the worst case, if they are stored. We have applied our algorithm on four 
kinship networks brazilian Indian ethnic groups: arara, deni, enawene- 
nawes and irantxe-myky. The results and consequences in Anthropology 
will be analyzed in a forthcoming paper. 
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Appendix 

We present here a proof of the Lemma [6j 

Lemma 6. Suppose a call to Single-Junction- All-Pairs(s, Sj), where 
Si € 6^ [s] . In line 3, the following assertions are true for the parameter 
vertex z and for all pairs of vertices u,v G T^. , such that post[u] , post[v] > 
post[w] + 1. 

— The vertex z is a representative; 

— The values p[u] and p[v] point to a representative; 

— For all vertices x £ Ts^\Tz, p[z] ^ p[x]; and 

— The vertex s G J{u,v) if, and only if, p[u] ^ p[v]. 

Proof. The first assertion is true as in the first call z = Sj is a representa- 
tive and before any recursive call to Single- Junction- All-Pairs(s, w) 
we make p[w] = w. The second assertion is true as line 5 or line 8 is 
executed on all w proper descendants of Sj and, in both cases, p[w] points 
to a representative. The third assertion is true as for any vertex v G Ts, 
p[v] points to a representative that is its ancestor in Tg. By first assertion, 
p[z] = z. But the vertex z is not an ancestor of x. Therefore, p[z] ^ p[x]. 
We prove the fourth assertion by induction in the number of iterations 
of the line 3. In the first iteration of line 3 we have z = Si and w is 
the child of the greatest post-order value, post[w\ = post[si] — 1. This 
means that Si is the single vertex in Ts^ with post-order value greater 
than post{w\. By definition, s ^ J{si,Si) and p[si] = p[si]. Now suppose 
that the fourth assertion is true in the i-th iteration of line 3 of a call 
to Single- Junction- All-Pairs(s, z). Let us prove that the assertion is 
true in the (i + l)-th iteration of this same call. Consider the directed 
path in Tz from z to w, Z = {z = wq, . . . ,Wk-i = w} in the i-th. it- 
eration. This path was configured during previous iterations and, for its 
construction, we know that s ^ J{z, wi), I = 0, . . . , k — 2. By Property [U 
post[z],post[t] > post[w] + 1, for z and for all t € 6q{w). Fourth assertion 
is valid for pair z,t, so s € J{z,t) if, and only if, p[z] p[t]. We can apply 
Lemma [3] and we will have two possibilities. 

Case 1. The vertex s ^ J{z,t), for all t G Sq{w). So, p[z] = p[t], for 
all t G Sq{ui) and by Lemma EJ s ^ J{z,w). The line 8 is executed and 
we have p[w] = z. Take x € Tg. with post[x] > post[w] + 1 (taking a 
vertex x with post[x] = post[w] we have, by definition, s ^ J{w,w) and, 
of course, p[w] = p[w]). Consider y = LCAt^{uj,x) (possibly x = y). 
Let us show that exactly before line 9 be executed, s E J{x, w) if, and 
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only if, p[x] / p[w]. Suppose x E Ts^\Tz. Then, the vertex y is some 
vertex of the directed path from Sj to z in T^. and, by third assertion, 
p[y\ 7^ p[A ^-iid p[x\ ^ p[z\ = z = p[w\. Therefore, we have to show that 
s G J(x, w). We have s € J(y, z) since / ^[z]. As G is a DAG, we can 
use the directed path from z to ii; in T^. to obtain s G J{y,w). If y = x, 
then is done. Else, by Lemma [H s G J{x,w). Now suppose x G Tj. So, 
y = LC^Tslw,^;) = Wj, wj G Z (possibly x = y). 

Case 1.1. The vertex s ^ J{x,z). In this case, p[x] = p[z]. Moreover, 
p[w] = z = p[z] = p[x]. So, we have to show that s ^ J{x,w). As we 
have said, s ^ J{z, wi), for / = 0, . . . , /c — 1. By Lemma El in particular 
s ^ J{y,w). Thus, if X = y, then s ^ J{x,w). li x ^ y, then, by the 
execution of the algorithm and Lemma [5l s cannot be a junction of any 
pairs X, Xr, where vertex in the directed path from z to x in T^. 

One of these vertices is y. Therefore, s ^ J{y,x). Finally, by Lemma 2, 
s ^ J{x,w) because s ^ J{y,w) and s ^ J{y,x). 

Case 1.2. The vertex s G J{x,z). Then, p[x] ^ p\z\. Moreover, p\ui\ = 
z = p[z] / p[x]. So, we have to show that s G J{x,w). We know that 
s ^ J{z,y). Then our assertion says ^[z] = So, p[y] = p[z] ^ p\x\. 
It means that s G J(y,x). Again, by Lemma 2, s G Jix^vS). After the 
execution of line 9, our assertion is restored. 

Case 2. The vertex s G J{z,t), for some t G Sq^w). So, 7^ for 
some 4 G Sq{w) and by LemmaO s G J{z,w). Line 5 is executed and we 
have p\w\ = w. Take x G T^. with post[x] > post[w\ + 1 (again, taking 
a vertex x with post[x] = post[vo\ we have, by definition, s ^ J{w,w) 
and, of course, = We are going to show that exactly before 

the execution of line 6, s G J{x,w) if, and only if p[x] ^ p[w]. As w 
is not an ancestor in Tg of x and p[x] points to a representative that 
is its ancestor in T^, we have p[x] 7^ p[w]. Thus, we have to show that 
s G J{x,w). If X G T2 \ Tw, then Lemma [Hb ensures s G J{x,w). If 
x £ Ts^\Tz, then y = LCAt^{x, w) is a vertex in the directed path from 
Si to z in Ts- (possibly x = y). By third assertion, p[y] ^ p[z]. Thus, 
s G J{y,z), and we can use the directed path from 2; to tt; in Tg- to 
obtain s G J{y,w). If x = y, then is done. If, x 7^ y, then, by Lemma 
El s G J{x,w). After the recursive call execution in line 6, we have for 
all G Tg., post[u],post[v] > post{YERTEx{minpost[w])) , the vertex 
s G J{u,v) if, and only if, p[u] ^ p[v]. Our assertion is restored after line 
7. □ 



